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Abstract — The constantly increasing demand for interactive 
broadband satellite communications is driving current research 
to explore novel system architectures that reuse frequency in a 
more aggressive manner. To this end, the topic of dual satellite 
systems, in which satellites share spatial (i.e. same coverage area) 
and spectral (i.e. full frequency reuse) degrees of freedom is 
introduced. In each multibeam satellite, multiuser interferences 
are mitigated by employing zero forcing precoding with realistic 
per antenna power constraints. However, the two sets of users that 
the transmitters are separately serving, interfere. The present 
contribution, proposes the partial cooperation, namely coordina- 
tion between the two coexisting transmitters in order to reduce 
interferences and enhance the performance of the whole system, 
while maintaining moderate system complexity. In this direction, 
a heuristic, iterative, low complexity algorithm that allocates 
users in the two interfering sets is proposed. This novel algorithm, 
improves the performance of each satellite and of the overall 
system, simultaneously. The first is achieved by maximizing the 
orthogonality between users allocated in the same set, hence 
optimizing the zero forcing performance, whilst the second by 
minimizing the level of interferences between the two sets. 
Simulation results show that the proposed method, compared to 
conventional techniques, significantly increases spectral efficiency. 



I. Introduction 

Towards the next generation of broadband multibeam Satel- 
lite Communication (SatCom) systems, innovative system ar- 
chitectures need to be considered in order to meet the highly 
increasing demand for throughput and close the digital divide. 
Spectrum scarcity is a major obstacle, especially in a satellite 
context where the higher frequency bands exhibit challenging 
channel impairments. In this direction, the investigation of 
full frequency reuse techniques that exploit the spatial degrees 
of freedom offered by the multibeam antenna is necessitated. 
Furthermore, in the evolution of geostationary (GEO) satellite 
systems, orbital slot congestion is an uprising problem. The 
fact that unpredictable changes in the traffic demand might 
cause the launch of secondary satellites that support existing 
ones should be considered as well. More importantly, long 
periods of coexisting satellites appear during the replacement 
phase of old satellites. Finally, cooperation between multibeam 
satellites has the potential of overcoming the major issue of 
Adjacent Satellite Interferences (ASI). It should also be noted 
that current research in cognitive satellite communications is 
considering scenarios where two satellites coexist [1]. All 



the above arguments provide reasonable cause towards the 
investigation of the optimal coexistence of two satellites. 

Existing literature on Multiuser Multiple Input Multiple 
Output (MU-MIMO) antennas, provides transmitter techniques 
to alleviate the multiuser interferences in many cases |2|. Fur- 
thermore, linear low complexity techniques have proven more 
realistic in terms of practical implementation. Specifically, 
focusing on the forward link (FL) of multibeam SatComs, 
linear joint processing techniques have shown great potential 
by providing a substantial tradeoff between implementation 
complexity and near optimal performance in terms of Sum- 
Rate (SR) [3|. Specifically, Zero Forcing (ZF) precoding, 
in detail analyzed further on, is based on channel inversion 
and then optimal power allocation over the users, with the 
aim of maximizing some performance metric. The metrics 
commonly addressed in literature involve either the total 
throughput performance (i.e. max SR criterion) or the Signal 
to Interference plus Noise Ratio (SINR) level of the worst 
user (i.e. max fairness criterion). Another important parameter 
in linear precoding is the type of constraints that will be 
assumed. Usually, a total sum power constraint simplifies the 
analysis and provides better results since the available power 
is freely allocated in every antenna. Despite its performance, 
the sum power constraint is unrealistic since each satellite 
antenna is fed by a dedicated high power amplifier (HPA) 
operating close to saturation. Hence, power cannot be allocated 
freely amongst the transmit antennas and a per antenna power 
constraint should be considered. In the present work, each 
satellite employs ZF beamforming while the power allocation 
is optimized under max SR criteria, subject to realistic per 
antenna constraints. 

Albeit the aforementioned technique, in a dual satellite 
system, intersatellite interferences still need to be handled. 
To completely mitigate these interferences, without frequency 
orthogonalization, full cooperation between the two trans- 
mitters has to be employed. Subsequently, the cooperating 
multi-antenna transmitters should perform joint or coherent 
transmission to all users while, both data and channel state 
information (CSI) must be exchanged flU. In a SatCom con- 
text, each satellite is served by one or multiple dedicated 
gateways, thus a fully cooperative dual satellite system would 
require a large number of interconnected GWs that exchange 



a substantial load of information. In the light of the above 
observations, partial cooperation (i.e. coordination) between 
the two satellites is proposed, therefore reducing the amount 
of data exchanged. 

The rest of the present paper is structured as follows. A 
brief review of the existing related work is provided in Section 
llll The considered system model is described with detail in 
Section Hn] Section ITVl explores the theoretical aspects of the 
developed heuristic algorithm, explaining the intuition behind 
it. In Section [V] simulation results for the performance of the 
proposed algorithm are presented. Conclusions are drawn in 
Section [Vll along with future extensions of the present work. 

II. Related Work 

Linear precoding techniques and especially ZF have been 
extensively investigated for terrestrial systems in Q, fl6j and 
the references therein. Despite their suboptimality, these tech- 
niques can still achieve asymptotically optimal performance 
under specific conditions, as proven in (|6)-||8). Specifically, if 
user channels are perfectly orthogonal to each other, ZF will 
attain maximum performance. Under the assumption of large 
random user sets, the probability of orthogonal users increases. 
User selection, however, is a highly complex problem. Despite 
its complexity, simple suboptimal algorithms in the existing 
literature provide substantial gains with affordable complexity. 
Based on existing algorithms, Yoo et al Q proposed an 
iterative user selection algorithm that allows ZF to achieve 
the performance of non-linear precoding j£| when the number 
of available for selection users grows to infinity. 

Despite the extensive literature on linear precoding and user 
selection, herein we further consider the optimal allocation 
of the selected users in two coexisting groups that interfere 
with each other. The novelty of the proposed work lies in 
the fact that user selection and allocation, not only optimizes 
the ZF performance of each system but also considers the 
interaction between the two transmitters, i.e. the inter-satellite 
interference. The procedure of selecting users out of a large 
pool and allocating them to specific sets is referred to as user 
scheduling. 

Notation: Throughout the paper, £[■], (•) , || • || denote the 
expectation, the conjugate transpose and the norm operations, 
respectively. Bold face lower case characters denote column 
vectors and upper case denote matrices while I„ denotes an 
identity matrix of size n. Upper case calligraphic characters 
denote sets. Operations A = , and A — B define an empty 
set and the relative complement of B in A, respectively. Finally 
\A\, denotes the cardinality of a set. 

III. System Model 

The system under investigation consists of two collocated 
multibeam satellites with overlapping coverage areas, serving 
fixed single antenna users (Fixed Satellite Services, FSS). A 
large number of users uniformly distributed in each beam 
is assumed and a Time Division Multiple Access (TDMA) 
scheme is realized, leading to one user per beam served, 
during each time slot. An overview of the considered system is 




Fig. 1. System model 



depicted in Fig.[T| where the focus is on the Forward Link (FL) 
downlink of the satellites (i.e. the link between the satellite 
and the users), while the FL uplink, or feeder link (i.e. the 
link between each satellite and the earth gateway station), is 
considered ideal. It should be clarified that the present work 
for simplicity purposes, does not model the earth curvature 
and the satellite orbit geometry. Subsequently, the variations 
in the distances of the beam centers and the distance between 
the satellites are not modeled but will be handled in future 
extensions. 

Considering each multibeam satellite separately, linear pre- 
coding is employed to cancel out multiuser interferences. By 
denoting N and K the number of transmit antennas and 
single antenna users respectively, the one user per beam per 
timeslot assumption implies a symmetric system, i.e. N = K. 
Subsequently, in each separate transmitter, a Multi-User (MU) 
Multiple Input Single Output (MISO) Broadcast Channel (BC) 
is realized and the input-output analytical expression for the 
fc-th user reads as 

y k = hj.x + n k , (1) 

where ht is an 1 x N vector composed of the channel 
coefficients between the A:-th user and the N antennas (i.e. 
feeds) of the satellite, x is an N x 1 vector of transmitted 
symbols and n k is the independent identically distributed (i.i.d) 
zero mean Additive White Gaussian Noise (AWGN) measured 
at the fc-th user's receive antenna. The noise is assumed 
normalized, thus £ {|nfe| 2 } = 1. 

To accurately model the multibeam satellite channel the 
following considerations are made. Under the assumption of 
fixed users with highly directive antennas, the satellite channel 
can be modeled as AWGN channel with real channel gain£] 
that depend only on the multibeam antenna pattern and on the 
user position. The elements of the fc-th user's channel vector 
are the square roots of the gain coefficients calculated using 
the well accepted method of Bessel functions iflOl 

la \ n ( Jl ( u ) i oa Mu) Y 

9ik(0ikj = &max — Ty 1" 36 — I ) (2) 

where u = 2.07123sin#/sin# 3( 2,B, J\, J 3 are the Bessel 
functions of the first kind, of order one and three respectively 
and G max is the maximum axis gain of each antenna. The 

'The main conclusions of this study, still apply if a random phase is 
incorporated in the channel gains. 



fc-th users' position corresponds to an off-axis angle On- with 
respect to the boresight of the i-th beam where Oi = 0°. 

A. Transmit beamforming 

Transmit beamforming is a multiuser precoding technique 
that separates user data streams in different beamforming 
directions [6j. Let us denote and p k as the unit power 

symbol, the N x 1 transmit precoding vector and the power 
scaling factor respectively, corresponding to the fc-th user. The 
total transmit signal will then read as 



K 



(3) 



k=l 



Subsequently, when beamforming is employed, (fTJ will be- 
come 



Vk 



hi 



1 ' 



rife, 



(4) 



where the first term of the summation refers to the useful 
signal and the second to the interferences. The column vector 
Wfe with 1 x N dimensions is the fc-th user's precoding vector, 
that is the fc-th column of a total precoding matrix W = 
[wi, W2, . . • , Wfc]. Subsequently, the SINR at each user is 



SINR,. = 



Pfc|hlw fc | 



1 



(5) 



When beamforming is employed, determining the optimal 
precoding vectors and power allocation vectors is tedious in 
practice. In the following paragraph a very common practical 
scheme, that will be employed in this paper is presented. 

B. Zero Forcing beamforming 

The implementation complexity of the MIMO BC channel 
capacity achieving dirty paper coding |9|, led to the devel- 
opment of less complex, yet suboptimal techniques. A linear 
precoding technique with reasonable computational complex- 
ity that still achieves full spatial multiplexing and multiuser 
diversity gains, is ZF precoding J6), ifTTl . lfl2l . The ability 
of ZF to fully cancel out multiuser interference makes it 
useful for the high SNR regime. However, it performs far 
from optimal in the noise limited regime. Also it can only 
simultaneously serve at most equal to the number of transmit 
antennas, single antenna users. A common solution for the ZF 
precoding matrix is the pseudoinverse of the K x N channel 
matrix H = [hi, h2 . . . , hj . Under a total power constraint, 
channel inversion is the optimal precoder choice in terms of 
maximum SR and maximum fairness [5|. However, according 
to the same authors, when Per Antenna power Constraints 
(PAC) are assumed, optimization over the parameters of a 
generalized inverse has to be performed. For simplicity, the 
present work only considers the simplified channel inversion 
design. Finally, due to the symmetrical system, the precoding 
matrix will read as W = H _1 . 

Following the complete mitigation of interferences, the 



SINR at the fc-th user will read as SINR 



ZF 



Pfc|htw fe 



Subsequently, the problem reduces to a power allocation 
optimization process, where the aim is to maximize the total 
system SR, subject to per antenna power constraints. This 
problem is formalized as 



max 

{p fc } 



s.t. 



K 



^log 2 (l + SINRf ), 

fc=i 

K 

fe=l 



(6) 



where Qjj is an K x K zero matrix, with the jj-th element 
equal to one and Pj = P to t/K, for every j, since we assume 
that all satellite antennas have identical RF chains. Problem 
©, is a standard convex optimization problem [13], |14| that 
is solved using common convex optimization tools 11151 . 

C. Dual satellite systems: channel model 

Extending the previous considerations to a dual satellite 
scenario, the resulting channel needs to be modeled. In this 
direction, two overlapping clusters of N\ and iV2 spot-beams 
covering K\ and K% fixed user terminals respectively, are 
considered (Fig. [TJ . The users, each equipped with a single 
antenna, are uniformly distributed over the coverage area. 
Despite the fact that in each satellite separately, a MU MISO 
BC is realized, the whole system operates over an interference 
channel. The fc-th user now has two vector channels, one 
towards each satellite, denoted using the indices 1 and 2 
respectively. The channel vectors hj^h]^ are the rows of 
a total channel matrix H tot of (K\ + K 2 ) x (JVi + N 2 ) 
dimensions, that models the satellite antenna gains of the two 
satellites. Assuming that the fc-th user is allocated to the first 
satellite its SINR will read as 



SINR,, = 



1 



EjL 2 iP2j|hi 2 w j2 | 



(7) 



In (0, the interferences from the adjacent satellite are clearly 
noted, while inter-satellite multiuser interferences are com- 
pletely mitigated by the the precoding. The equivalent relations 
for the users allocated to the second satellite, are straightfor- 
ward to deduce, by exchanging indices 1 and 2. 

IV. User Selection and Allocation 

As proven in Q, user selection can significantly improve 
the performance of ZF in an individual system. However, 
considering the coexistence of two separate transmitters, as 
is the case in a dual satellite system, partial cooperation, 
namely coordination, can be employed to solve the problem of 
high intersatellite interferences. To this end, this contribution 
proposes an algorithm that selects users and allocates them 
to each satellite. Intuitively, the two basic criteria that need 
to be considered for this procedure are: a) the maximization 
of the performance of each satellite separately, and b) the 
minimization of the interferences between the two sets. 

The performance of each satellite separately is optimized by 
constructing a semi-orthogonal user group from a vast number 
of users JS). Extending this result, the creation of two user 



sets under the semi-orthogonality criterion is straightforward 
since the channel gain of each user can be projected to the 
orthogonal complement of the channels of the previously 
selected users. In each iteration of the algorithm, the user with 
the maximum projection is allocated to the corresponding set. 
This simplistic approach has been considered for comparison 
reasons. 

The novel proposed algorithm accounts for the effects of the 
interferences between the two sets. It should first be mentioned 
that the exact calculation of the level of interferences in 
each iteration is not possible since the exact user set is 
still undetermined. To exactly calculate the interferences, one 
would need to solve the power optimization problem (0 for 
all possible combinations of users. Under the assumption 
of large number of users, this would lead to unaffordable 
computational complexity. However, based on a basic ad- 
vantage of ZF beamforming, which is the decoupled nature 
of the precoder design and the power allocation optimiza- 
tion problems, an approximation of the interferences can be 
made. In this direction, the precoding vectors of the users 
selected in the previous iterations can be utilized to provide 
an indicative measure of the interferences between the user 
sets. This implies that an equal power allocation is assumed. 
This assumption becomes asymptotically exact in the high 
SNR regime, where the powers allocated to each user are 
approximately equal. Incorporating all the above, a heuristic, 
iterative Semi-orthogonal Interference aware User Allocation 
algorithm (SIUA) has been developed and will be presented 
in the following. 

The SIUA algorithm, presented in full detail in Algo. Q] 
works as follows. During the initialization procedure, i.e. Step 
1, the strongest user towards each satellite is allocated to the 
equivalent group. While the two sets are not full, Steps 2 
and 3 are executed. In Step 2, for each of the unallocated 
users, the following metrics are calculated: a) In accordance 
to 0, gife and g2fc represent the orthogonal component of 
each unallocated user's channel to the orthogonal subspace of 
the already allocated users, for the two sets respectively. In b) 
IL and V k2 are equivalent measures for the interference each 
user would receive if equal power allocation is employed. It 
is calculated as the squared norm of the product of the users' 
channel with the power of the transmit signal of the second 
user set and the channel of each user. Finally, in c) V kl and 
P k2 are approximations of the interference that the allocation 
of each user can potentially induce to the second set, if this 
user is allocated in the respective set. It is calculated as the 
product of the interferences this user induces to every user that 
belongs to the second set. Since the goal is to find the most 
orthogonal users that at the same time receive and induce the 
least possible interferences, the measure to be maximized is 
the fraction of the orthogonality metric over the product of the 
interference metrics. At the last stage of each iteration, two 
maximum fractions /ii and n% are calculated over the hole user 
set and compared between them. The user that corresponds 
to the largest measure among the two, is allocated to the 
equivalent satellite. 



TABLE I 
Link Budget Parameters 



Parameter 


Value 


Orbit 


GEO 


17-,.,.. ■ ,.. . U„„ A 

frequency band 


Ka (20 GHz) 


User link bandwidth 


JUL) MriZ 


I'lLllllUCl Ul UCttlllS 


7 


Beam diameter 


600 Km 


TWTA RF power @ saturation 


+21 dBW 


Max satellite antenna gain Gt 


+52 dBi 


Max user antenna Gain Or 


+40 dBi 


Free space loss 


-210 dB 


Signal Power S 


-97 dBW 


Receiver noise power N 


-118 dBW 


SNR S/N 


21 dB 



The described heuristic, iterative, optimization algorithm, 
only requires full knowledge of the CSI of all users. Conse- 
quently, coordination reduces the amount of data that needs 
to be exchanged since each GW handles only the data of 
the users allocated to the corresponding satellite. Moreover, 
SIUA runs only as many times as the number of transmit 
antennas and thus compromises a scalable solution that can 
be extended for larger multibeam systems. Another advantage 
of this solution is that the power optimization in each satellite, 
a convex optimization problem that requires some computa- 
tional complexity, is decoupled from the algorithm execution. 
Additionally, despite the fact that the solution is heuristic and 
not optimal, it is considerably less complex since the optimal 
user allocation would require exhaustive search of all possible 
combinations of the users. Finally, SIUA, can be executed in 
a centralized location or run in parallel at the GWs that share 
CSI. 

V. Simulation results 

To evaluate the performance of the proposed algorithm, two 
satellites each with seven beams where assumed. The low 
number of beams is only chosen to reduce the simulation time 
of the convex optimization problem © and has no effect on 
the algorithm, as discussed in Sec. [IV] It is however inline 
with the future considerations for the terabit satellites, where 
each GW is expected to handle between 5 and 8 beams. 
Additionally, the simulations are performed according to the 
link budget calculations described in Table |IJ where it can be 
noted that the normal SNR operating point of current satellite 
systems is 21dB. 

In Fig. |2] the results of Monte Carlo simulations that 
calculate the performance in terms of system Sum Rate (SR) 
of the coexisting systems with and without cooperation, are 
presented. Due to random user positioning, 100 iterations 
where executed, each with a different user position pattern 
so that the average performance can be evaluated. The upper 
bound for the system performance is deduced assuming full 
cooperation amongst the transmitters. For this case two curves 
show the SR in bps/Hz: one for the average performance of 
random user positioning and one employing the algorithm 
developed in J6|, namely the Semi-orthogonal User allocation 
(SUS) algorithm, which allocates users without regarding the 



coexistence of the systems. Subsequently, an average gain of 
25% is noted by employing a simple user selection scheme, 
instead of randomly selecting users. Since in the fully coop- 
erative system, interferences are completely mitigated, SIUA 
is unnecessary. 

The substantial performance gain from SIUA is proven in 
Fig.[2]for the more realistic case of coordinated systems. In this 
figure, SIUA is compared to SUS and also to an independent 
interfering system. From these curves, it is concluded that in 
the low SNR region, the SUS algorithm performs better, as 
expected since the noise limited regime is almost interference 
free. In the SNR region of interest, however, it is clear that 
the SIUA algorithm, by reducing the level of interferences, 
provides substantial gains: more than 52% of improvement in 
terms of SR, over a non cooperative system and 28% of im- 
provement over a coordinated system employing simple user 
selection. It is therefore concluded that the SIUA sacrifices 
some low SNR performance to provide substantial gains in the 
SNR region of interest. A simple switching scheme between 
the two algorithms can provide good performance over the 
whole SNR region of interest. 

In Fig. [3] a coordinated system using SIUA, is compared 
to an ideal non interfering dual satellite system that em- 
ploys frequency orthogonalization to allow the operability 
of the two coexisting satellites. This approach models the 
currently employed techniques of bandwidth splitting. In this 
plot it is proven that around the SNR area of interest, the 
proposed algorithm outperforms the conventional techniques 
(25% gain). Therefore, the SIUA comprises a candidate tool 
for handling ASI. As the SNR increases, the gain decreases 
as expected, since the conventional system operates under the 
ideal assumption of zero interferences. 

Finally, in Fig. [4] the behavior of the algorithm with respect 
to the size of the user pool is investigated. To this end, 
the achievable SR for a given value of SNR, i.e. 20dB, 
is calculated as the total number users increases and also 
compared to the performance of the SUS algorithm. In this 
figure it is proven that the algorithm reaches close to its 
maximum performance for a finite number of total users (600 
users) and further increase of the user pool has little effects. 
From the same figure we note that the rate of convergence of 
the proposed technique to the saturation point is very similar 
to the SUS algorithm. 

VI. Conclusions and future work 

The topic of dual satellite systems is addressed in the present 
contribution and a suboptimal, simple solution is proposed. 
A low complexity, heuristic algorithm that minimizes the 
interferences between the two groups, while maintaining the 
orthogonality between the users of the same group, allows 
for the coexistence of two separate multibeam, joint pro- 
cessing, coordinated satellites. Thereby the overall system 
spectral efficiency is increased. The only cost of the proposed 
solution is the CSI exchange between the GWs. According 
to simulation results, the algorithm achieves 52% of spectral 
efficiency improvement over non-cooperative full frequency 
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Fig. 2. Evaluation of SIUA algorithm in terms of system sum rate, by 
comparison with optimal and interfering systems. 

reuse systems and 25% improvement over non-cooperative 
conventional orthogonalized in the frequency domain systems. 
Subsequently, the proposed scheme, successfully exploits the 
spatial orthogonalization of users and allows for systems that 
are partially cooperative to operate over all the available 
spectrum. 

This fundamental approach can be extended as a useful 
interference mitigation technique for several other scenarios. 
Also, the consideration of a third set of users that will 
be simultaneously served by both satellites to achieve some 
minimum requirements will be part of future work. The dual 
satellite system will also be investigated as a cognitive com- 
munications system where priority will be given to primary 
users. 
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